EPJ manuscript No. 

(will be inserted by the editor) 



Thermalisation by a boson bath in a pure state 

S. Camalet 

Laboratoire de Physique Theorique de la Matiere Condensee, UMR 7600, Universite Pierre et Marie Curie, Jussieu, Paris-75005, 
France 

Received: date / Revised version: date 

Abstract. Wc consider a quantum system weakly coupled to a large heat bath of harmonic oscillators. It 
is well known that such a boson bath initially at thermal equilibrium thermalises the system. We show 
that assuming a priori an equilibrium state is not necessary to obtain the thermalisation of the system. We 
determine the complete Schrodinger time evolution of the subsystem of interest for an initial pure product 
state of the composite system consisting of the considered system and the bath. We find that the system 
relaxes into canonical equilibrium for almost all initial pure bath states of macroscopically well-defined 
energy. The temperature of the system asymptotic thermal state is determined by the energy of the initial 
bath state as the corresponding microcanonical temperature. Moreover, the time evolution of the system 
is identical to the one obtained assuming that the boson bath is initially at thermal equilibrium at this 
temperature. A significant part of our approach is applicable to other baths and we identify the bath 
features which are requisite for the thermalisation. 

PACS. 05.70.Ln Nonequilibrium and irreversible thermodynamics - 05.30.-d Quantum statistical mechan- 
ics 



1 Introduction. 

The environment of a physical system plays an important 
role in its dynamics. The environmental degrees of free- 
dom are responsible for the relaxation of the system into 
thermal equilibrium. As is well known, this irreversible 
evolution follows from Schrodinger dynamics for a sys- 
tem weakly coupled to a large heat bath, provided the 
bath is supposed to be initially at thermal equilibrium. 
However, this way to model the thermalisation process 
leaves open the question of how the system environment 
reaches its own equilibrium state. An equilibrium state 
is also assumed a priori in the most common derivation 
of the canonical equilibrium distribution. In this standard 
argumentation, the interaction between the system and its 
heat bath is neglected jl] and the microcanonical princi- 
ple [5] is invoked, i.e., the composite system consisting of 
the considered system and the bath is assumed to be in 
a microcanonical mixed state. The ensuing reduced state 
of the subsystem S of interest is a canonical thermal mix- 
ture of its self-Hamiltonian eigenstates. Recently, it has 
been shown [3,4.5 that no statistical averaging is actu- 
ally required to derive the canonical distribution for S. 
This thermal state is found for almost all pure states of 
the composite system of macroscopically well-defined en- 
ergy. Canonical equilibrium is thus obtained as a direct 
consequence of the quantum mechanics principles. 

Nonetheless, to understand the thermalisation process, 
it is necessary to take into account the interaction between 
S and the bath and to study the resulting Schrodinger dy- 



namics of the isolated composite system. For the special 
case where the system-bath coupling Hamiltonian com- 
mutes with the self-Hamiltonian Hs of S, it has been 
shown [6l[7] that S evolves into a statistical mixture of 
the eigenstates of Hs for an initial "generic" pure product 
state of the composite system. But such a coupling cannot 
lead to a complete thermalisation of S. As the Hamilto- 
nian Hs is a constant of motion, the populations of the 
eigenstates of Hs remain equal to their initial values. More 
recently, some results have been obtained regarding the 
thermalisation induced by a heat bath initially in a pure 
state for abstract system-bath interaction Hamiltonians 
defined in the eigenbasis of the non-interacting Hamilto- 
nian. Analytical arguments showing that S is at equilib- 
rium most of the time are presented in [3j and examples 
are studied numerically in [5]. 

In this paper, we consider a system S weakly coupled 
to a heat bath of harmonic oscillators. Numerous physi- 
cal environments are well modelled by such a boson bath 
which is thus often used to describe dissipation in quan- 
tum systems [9] . From this perspective, the bath is usually 
assumed at thermal equilibrium at initial time. Here we 
study the reduced dynamics of S for an initial pure prod- 
uct state of the composite system consisting of S and its 
bosonic environment. Using the resolvent operator formal- 
ism, we derive this dynamics from the Schrodinger equa- 
tion of the isolated composite system. Since we are inter- 
ested in the thermalisation process, we suppose S is ini- 
tially in a pure state. For the bath, we consider a "typical" 
state of macroscopically well-defined energy. Whereas our 
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results are obtained for a boson bath, a significant part 
of the formalism developed is applicable to other environ- 
ments. The paper is organised as follows. In the following 
section, we present the model and the initial bath state. 
The reduced dynamics of the subsystem S is derived in 
Sect. [3] We then determine the asymptotic state of S and 
the relaxation to this state in the weak coupling regime in 
Sect. The results and their potential validity for other 
environments are discussed in Sect. [3 

2 Model 

We write the total Hamiltonian of the composite system 
consisting of the system S of interest and its bosonic en- 
vironment as 

H = Hs + Hb+Hj (1) 

where H$ and Hi describe, respectively, the intrinsic dy- 
namics of S and its coupling to the boson bath. The di- 
agonalised bath Hamiltonian H E reads 

H B = J2^ b lK ( 2 ) 

9 

where b q (b q ) are bosonic creation (annihilation) operators 
and the sum runs over the N harmonic modes of the bath. 
The eigenfrequencies w q depend on the specific bath con- 
sidered [jJJ. The thermodynamic limit TV 3> 1 is assumed 
throughout this paper. In the following, the eigenstates 
of Hs and Hb are denoted by Latin and Greek letters, 
respectively, i.e., 

H s \k)=e k \k) 

H B \a) = E a \a) (3) 

where e& and E a are the respective eigenvalues. For the 
system S, we consider a discrete spectrum of nondegen- 
erate eigenergies efc. A system with continuous spectrum 
cannot relax into thermal equilibrium. For example, in the 
case of a free damped particle, the width of an initially lo- 
calised state grows indefinitely. The bath influences signif- 
icantly this diffusion process [10] . For the bath, an eigen- 
state | a) corresponds to a set of harmonic mode occupa- 
tion numbers {n q } and E a — ^2 q <-<J q n q . In the thermo- 
dynamic limit, the eigenenergy spectrum of the bath can 
be characterised by a density of states n(E) given by the 
standard following definition [5] . For a macroscopic energy 
E ~ N and an energy SE <C E but far larger than the 
maximum level spacing of Hb, the number D of states \a) 
such that E < E a < E + SE, is practically proportionnal 
to SE and n(E) = D/5E. Moreover, this density satisfies 
the Boltzmann's relation 

In [n(E)SE] ~ Ns(E/N) (4) 

where s is the bath entropy per oscillator. We use the units 
h = Ub = 1 throughout this paper. As the bath entropy 
Ns and energy E are extensive variables, the entropy per 
oscillator s assumes finite values for finite energies per 
oscillator E/N. 



The interaction Hamiltonian Hi depends on the de- 
grees of freedom of both the bath and the system S. It is 
useful for the following to expand it as 

^ =j2m\u k i (5) 

where Uu are bath operators which obey lit, = U lk . This 
form is generally valid for a composite system. We consider 
here the interaction operators 

Ukl =J2<( b l + b q )- (6) 
i 

Numerous environments can be modeled by such a linear 
coupling to the positions of independent harmonic oscil- 
lators [9]. The influence on S of the bosonic environment 
described by Hb and the coupling operators ((6]) can be 
characterised by the standard spectral densities 

Jij,*l(«) = K qffi(w - U q ) (7) 

9 

which are finite continuous functions of u> > in the ther- 
modynamic limit. Since (n 3 q )* — k 1 ^ , they satisfy Jkl,ij 

J ■ T 

J ji,lk ~ J ij,kl- 

In this paper, we study the Schrodinger time evolution 
of the system S following from an initial pure product 
state of the composite system. We assume that the initial 
bath state 

iv) = y, ^ a > ( 8 ) 

is a normalised state of the Hilbert space H e defined by 
the D states \a) such that E < E a < E + SE. In the 
following, we show that almost all normalised states in 
He lead to the same behavior of S. To do so, we use the 
uniform measure \i on the unit sphere in He given by 

A*({iM = ^zr^(i- E l^l 2 )- ( 9 ) 

\a)en E 

With this normalised measure, we will prove that a subset 
of bath states (jSJ) resulting in the same time evolution of S, 
is of size unity in the thermodynamic limit. In [5], the uni- 
form measure is discussed for composite system states <P 
of macroscopically well-defined energy. The product states 
considered here and the corresponding measure © can 
be obtained by performing ideal quantum measurements 
of the observable Hs on the states <P. We remark that 
other choices are possible to describe an initial bath state 
of macroscopically well-defined energy. For example, the 
space He can be replaced in © by the Hilbert space H' E 
defined by the eigenstates \a) such that E a < E. In this 
case, the appropriate measure is given by the expression 
© with the replacements He — ► H' E and D — » D' where 

D' = n{E')dE' is the dimension of H' E . Due to the 
exponential A^-dependence of the density of states n, the 
number of states \a) in H' E but not in He-se is negligible 
and hence the energy of a typical state in H' E is practi- 
cally equal to E ~ N. For the same reason, the results 
presented in the remainder of the paper are valid for such 
a state. 
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3 Reduced dynamics 

We derive here the reduced dynamics of the subsystem S 
for an initial disentagled state 



M = PSPB 



(10) 



where ps and ps are density matrices of S and the bath, 
respectively. The initial bath state can be, for example, a 
thermal mixed state, ps oc exp(— Hb/T), or a pure state, 
Pb = In this paper, we are concerned with the 

thermalisation induced by the pure states given by (JSj> . 
However, the calculation presented in this section is not 
restricted to pure initial bath states and it is interesting 
for the following to consider also thermal states. To obtain 
the time evolution of S, we first write its reduced density 
matrix at time t > as 



P(t) 



2vr 



dze~ lzt p(z) 



(11) 



where r\ is a positive real number. The Laplace transform 
p(z) of the state p(t) is given by 



p(z) = Tr s [Q(z)0] 



(12) 



where Tr^ denotes the partial trace over the bath degrees 
of freedom and Q{z) = (z — C) 1 is the resolvent of the 
Liouvillian C defined as CA = [H, A] for any operator 
A. The Liouvillians Cs, Cb and Ci corresponding to the 
Hamiltonians Hs, Hb and Hj, respectively, are given by 
similar definitions. An equation for p(z) can be derived 
from the Schrodinger time evolution of the composite sys- 
tem consisting of S and its environment using the projec- 
tion superoperators V and Q defined by Q = 1 — V and 
VA = Tr b (A) ps where A is any operator. With these 
superoperators, we deduce from (z — C)Q(z)Q — Q the 
coupled equations pTllT2] 

v(z - c)vvQ{z)n + v(z - c)QQQ(z)n = n 

Q(z - C)VVQ{z)Q + Q{z - C)QQS(z)n = (13) 

for the operators VQ[z)Q and QQ{z)£l. Solving the latter 
for QQ(z)Q in terms of VQ(z)f2 and then substituting in 
the former yields 

Tr B [(z-C- CQ(z - QCQ)~ l QC) p B p{z)] = Ps . (14) 

This equation describes the complete time evolution of 
S. Since we are interested in the limit of weak coupling 
between S and the bath, we expand (z — Q£<2) _1 in the 
Liouvillian Ci- With this expansion and the properties 
Ty b {£bA) = and Ti B {C s A) = £ s Tr B (A) where A is 
any operator, we find 



Ps 



zp(z) 
-Tr B 



[Hs,p(z)} - Tt b [C lPB p{z)] 
1 



(15) 



where Csb 



n>0 

Cs + Cb- 



C 



SB 



Q) (C B +Ci)p B p{z) 



It is now convenient to expand the Laplace transform 
p(z) in the basis \k)(l\ of the Liouville space of S as 



p(z) = ^2\k){l\r kl (z) 



(16) 



k.l 



where rki(z) are scalar functions of z. These components 
are solutions of the coupled equations 

(z - ^) rij {z) - E$(z)r M (z) = (i\p s \j) (17) 
k,l 

where u>ij = £j — tj. The influence of the environment is 
described by the self-energies 



E*{z) = <i|T ri 



X p B \k)(l\ 

z - C S b 11 



(18) 



To derive this expression from (fT5")) we have grouped the 
terms of same order in Hj using 1 + B^ 1 (A — B) =B~ X A 
with A = z — Cs and B — z — Csb- For the usually con- 
sidered thermal bath state ps oc exp(— Hb/T), the above 
derivation is simpler as CbPb = and hence Q and Csb 
commute. In this case, the superoperator (z — uiki)/{z — 
Csb) hi (f!8)) can be replaced by 1 [13] . We remark that, 
as Tr(£/^4) = for any operator A, the self-energies (fT8| 
satisfy the relations J2i ^u( z ) — which ensure that the 
trace of the density matrix p(t) is conserved. To obtain 
the reduced dynamics (fiT]) we have only used the form 
(HJ) of the total Hamiltonian describing the system S, the 
bath and their coupling. The bosonic nature of the bath 
and the specific coupling bath operators ((6|) do not play 
any role in the above calculation. 

The equations (|17p can be solved in rki(z) by writ- 
ing these components as power series in Hi. However, 
the perturbative solutions obtained by this method give 
the state p(t) of S, through (fT6| and (fTT|) . only for short 
times t. To determine the long time evolution of S for 
weak coupling to the bath, we proceed as follows. The 
coupled equations (fT7|) can be written as M(z)r(z) = rg 
where r(z) and rg are vectors of components r k i(z) and 
(k\ps\l), respectively, and M(z) is a matrix of elements 
Mij^i(z) = (z — ujij)6kiSij — U^J(z). This matrix equation 
can be formally solved by inverting M(z) with the help 
of its eigenvalues X n (z) and eigenvectors u n (z) and the 
eigenvectors v n {z) of the transpose M(z) T . With these 
definitions, we write the solution of (fT7|) as 



r(z) = 



1 



^ Xn(z) 



u n (z)v n (z) T r S - 



(19) 



In absence of coupling to the bath, i.e., for Hi = 0, the 
matrix M(z) is diagonal and {A n (z)} = {z — LUki}- As a 
result, the populations (k\p(t)\k) are constant and the co- 
herence (k\p(t)\l) oscillates with the frequency LOki- For a 
weak coupling, the eigenvalues A„(z) can be determined 
as power series expansions in Hi as discussed in the next 



4 



S. Camalet: Thermalisation by a boson bath in a pure state 



section. As a consequence of the property J2 i S^{z) = 0, 
Ao(^) = z remains an eigenvalue of M(z) for any coupling 
strength. Moreover, the components of the corresponding 
eigenvector v (z) are dij for any z. The term n = of (fl"9l) 
thus simplifies to Uq(z)/z. This term results in a steady 
contribution to the state p(t) of S. For the model consid- 
ered in this paper, any eigenvalue of index n ^ manifests 
a branch cut on the real axis and A n (w + i0 + ) is a non- 
vanishing function of the real argument oj. Consequently, 
the terms n ^ of (fl~9|) correspond to contributions to 
p(t) which vanish at asymptotic times. In conclusion, the 
components of the asymptotic state poo of S are that 
of u (i0 + ), i.e., they are given by 



k.l 



(20) 



Using this equation and the expressions (|18p . p^ can be 
determined for weak coupling to the bath by writing the 
components as power series in Hj. 



4 Weak coupling regime 

In this section, we determine the asymptotic state poo 
of the system S for weak coupling to the boson bath 
and a typical initial bath state |[8]) of macroscopic en- 
ergy E. We also discuss the terms n ^ of the expres- 
sion ([TO]) which describe the decoherence and relaxation 
of S towards p^. More precisely, we evaluate the low- 
est relevant orders in if / of the state p^ and eigenvalues 
X n (z). To do so, we first derive tractable power series ex- 
pansions of the self-energies (fl"8|) as follows. By writing 
(z — Csb)~ 1 as the Laplace transform of exp(— UjCsb) 
and Ci in terms of the bath operators U k i using ([5]), 
we express SfJ(z) in terms of the Heisenberg operators 
Uki(t) = exp(itHB)Uki exp(— itHs)- We will see that, for 
weak coupling to the bath, the reduced dynamics of the 
subsystem S is essentially determined by the averages 
(U k i(t)) and (Uij(t)Ukt(t + T)) where (A) = (V>|A|V) for 
any bath operator A. As shown below, explicit expressions 
can be obtained for these expectation values. 



4.1 Correlations of the bath interaction operators 

We show here that, in the thermodynamic limit, almost 
all bath states ((8|) lead to the same averages (Uki(t)) and 
(Uij{t)Uki{t + t)). We will obtain this result by comput- 
ing the Hilbert space averages [4] and variances of these 
expectation values resulting from the measure ([9|) . For ex- 
ample, the Hilbert space average of (Uki(t)) reads 

(M*I) = I1 / ^Vv»MWa})X>^/j(a|tf«(*)l/9) 

a J a ,/3 

(21) 

where d 2 ip a — dReip a dlmip a and the product and the 
sums run over the bath eigenstates \a) £ He- To evaluate 



this integral, we first note that, for a bath size N 3> 1, 

J] [d 2 ^({M)^f[^e- D M (22) 

ct^at(p) p=l 

where the n components ip a ( p ) are fixed [14] . This reduced 
distribution and the expansion ((5J) then yield (Uki(t)) = 
{Um)e where 



(A) E = ^ £ (a\A\a) 



(23) 



denotes the microcanonical average of any bath operator 
A. We define the variance a\ as the Hilbert space average 
of \(U kl (t))-(U kl ) E \ 2 and find 

^ E IH^I/5)| 2 <^(C/ fei ^> £ . (24) 

The upper bound is simply obtained by replacing the sum 
over the states \(3) £ He by a sum over all the bath eigen- 
states |/3). For (Uij(t)Uki(t + T)), we find the Hilbert space 
average (UijUkiir)) e and variance 



4(r) 



1 

£2 



E IH^,t4*(T)|/3}| 2 

|a),l/3)eW E 



(25) 



We now determine the microcanonical correlation func- 
tions (UijUkiir)) e using the equality of the occupation 
number (b q b q )E to the canonical average (b q b q )x with the 
appropriate temperature T. The canonical average of any 
bath operator A for temperature T reads 

1 V 



(Ah 



-E a /T 



a\A\a) 



(26) 



where Z = J2 a ex P( — E a /T). It can be established, by the 
following standard arguments [2], that (b q b q )E = (b q b q )T 
for T equals to the bath microcanonical temperature cor- 
responding to the energy E, i.e., 



T 



(27) 



where s' is the derivative of the bath entropy per oscilla- 
tor s with respect to the energy E/N. The microcanon- 
ical distribution of the harmonic mode occupation n q is 
P(n q ) = n(E—n q u} g )/n(E) where n is the density of states 
of the system consisting of the N — 1 other bath modes. 
This density satisfies the Boltzmann's relation (j4|) with 
the corresponding entropy s. Expanding this entropy in 
the energy n q oj q and taking into account that s — s in the 
thermodynamic limit, results in P{n q ) oc exp(~n q u) q /T) 
where T is given by (I27|) . Using this distribution, we ob- 
tain 



(UijU k i(T)) E = (UijUkiir))! 



du>Jij,ki(u) 



cos(cjt) 
tanh(w/2T) 



i sin(o;T) 



(28) 
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It follows from this result that <j\ = 0(exp(— N)). 
Consequently, the overwhelming majority of bath states 
© satisfy 

(U kl (t)) = 0. (29) 

The dynamical fluctuations of (Uki(t)) around zero disap- 
pear in the thermodynamic limit for typical bath states 
\ip). By generalising the above derivation to the micro- 
canonical distribution of two harmonic mode occupa- 
tion numbers, it can be shown that the left side of the 
inequality (|2"5|) is equal to a canonical average. More- 
over, as the operators U k i are linear combinations of 
the operators b q and b q , the upper bound of cr|(r) can 

be rewritten as ({U i3 Uji)T(UkiUik)T + \{UijUki{T))r\ 2 + 
\{UijUik{r))T\ 2 ) I D '. As a result, the variance of (r) also 
vanishes exponentially in the limit N 3> 1 and hence, us- 
ing (HU, 



(U t3 (t)U kl (t + T)) = (UijUuir)) 



T- 



(30) 



4.2 Asymptotic state, decoherence and relaxation 



We deduce from the result (f2"9")l that the first order con- 
tribution in Hj to the self-energy (fl8|) with ps = \ip) 
vanishes. To determine the second order contribution, we 
first rewrite it as a sum of four terms using QA = A — 
PbTtb(A) where A is any operator. As the first-order 
self-energies vanish and Tr^[(z — Csb) -1 PbTib(A)] = 
Ttb[{z — Csb) 1 A] where A is any operator, only one 
term remains and 



E$(z) = (z - ujki)(i\Tr B 



{c. 



1 



£-SB 



PB\k)(l\ 



(31) 

to second order in Hi, where ps = A lengthy but 

straightforward calculation then gives 



+5ji)^e lt ^(U ip (t)U pk }T 
p 

{Ui 3 {t)U lk ) T - e lt ^ k {U l3 U t k(t)h 



(32) 



to second order in Hi. The final expression of the calcula- 
tion has been rewritten with the help of the equality (|5U|) . 
The temperature T, given by (|27j) . is the microcanonical 
temperature of the bath corresponding to the energy E of 
its initial state. 

Since u>ij ^ for i ^ j and the self-energies Efj are at 
least of second order, (j2"0|) leads to — PkSki + O(Hj). 
Writing this equation to second order for i — j and using 
(l32l). we find 



^ r lk p k -YihP, o. 



fe 



(33) 



The zeroth-order populations Pi are thus the steady solu- 
tions of a master equation. From (|32|) and (|28|) . we deduce 



the (positive) rates 



Pik — 



K lk (e % — efc) 

e (ei-e k )/T _ l 



(34) 



where K ik (u) = 2n[0(v)J i i Sl ki(u) — G(—tjj)Jik,ki(—v)\. As 
Jki,ik = Jik,H> the rates (|34|) satisfy the detailed balance 
relation, 



P k i 6 ^ Pq 



fee 



-e k /T 



(35) 



Consequently, assuming the matrix G of elements Gi k = 
—Pik + SikPki is irreducible [2], the zeroth-order asymp- 
totic state of S is the canonical thermal state, 



-H S /T 



Tr(e- H s/T) ■ 



(36) 



The following special case is worth mentioning. If the self- 
Hamiltonian H$ of S and the interaction Hamiltonian Hi 
commute, i.e., U k i = for k ^ I, the populations (k\p(t)\k) 
remain equal to their respective initial values. In this case, 
the only effect of the environment is to destroy the coher- 
ences (k\ P {t)\i) nm. 

We now discuss the terms n ^ of the expression 
(jTOjl which describe the decoherence and relaxation of S 
towards the thermal equilibrium state (f36|) . We first con- 
sider a frequency Wki 7^ such that u>ij ^ Uki for any 
(hj) (^)O- F° r such a frequency, we find X n (z) = 
z — ujki — Euj (z) to second order in Hi. The corresponding 
term in (fl9|) describes the time evolution of the coherence 
(k\p(t)\l). The influence of the environment on this coher- 
ence is easier to understand for z = lu + i0 + where lu is 
real [TTlfl2] . For weak coupling to the bath, the perturba- 
tive corrections to \ n (u + iO + ) can be neglected except for 
to ~ Lo k i ■ From ([32]) and (|28|) , we obtain 



+ 



2^ 

p 



du>- 



2 ^ 

p^k,l 
1 

»"/ T - 1 



{P pk + Ppi) — iS, 



ki 



K pk (uj) K pl ( 



- u pk 



w p i 



(37) 

to second order in Hi, where p denotes the Cauchy prin- 
cipal value and S k i = ^J2 q ( K q k ~ K q) 2 coth(ui q /2T)S(oj q ). 
The negative imaginary part, in the first line, leads to an 
exponential decay of (k\p(t)\l) at long times (9l im[T2"] . The 
real part results in a slight shift of the oscillation frequency 
of this coherence. As is well known, there exists a short- 
time regime which is not accurately described by the ap- 
proximation X n (u! + iO + ) ~ lu — uiki — E^iuiki +i0 + ). How- 
ever, the behavior of S at short times t is also determined 
by the second-order self-energies (1521 as p(t) is well ap- 
proximated by an expansion to second order derived from 
(fT7|) . If LOki = (I" 7^ for several (k,l), the second-order 
contributions to the eigenvalues A„(z) = z—lu+- ■ ■ are ob- 
tained by diagonalising a matrix of elements E^j(0j + i0 + ) 
(to second order) where (i, j) and (fc, I) are such that coy = 
LOki = The corresponding coherences of p(t) vanish at 
asymptotic times. For the eigenvalues A„(z) = z + ■ ■ ■ , it 
can be shown that A„(iO + ) = i"/ n to second order, where 
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7„ are the positive real eigenvalues of the matrix G. One of 
them vanishes and corresponds to the asymptotic state of 
S given by ([33]) . The other ones are the characteristic rates 
of the long-time relaxation of the populations (k\p(t)\k). 
In conclusion, for weak coupling to the bath, the complete 
time evolution of S is determined by the second-order self- 
energies (|32|) and hence is identical to the one obtained in 
the usual case of a bath initially at thermal equilibrium, 
PB = exp(-H B /T)/Z. 



5 Discussion of the results 

We discuss here the relative importances of the different 
environment features in the obtention of the asymptotic 
thermal state (ffi)]) . We first reiterate that the derivation 
of the reduced dynamics of the subsystem S 1 , presented in 
Sect. [3l relies essentially on the form ([1]) of the total Hamil- 
tonian of the composite system. To obtain then from l|20[) 
the equations (|33[) for the zeroth-order asymptotic pop- 
ulations, the vanishing of the first order contributions to 
the self-energies (|18l) is not necessary. It can be proved 
from the expression (|18p that the equality (|31[) is correct 
to second order in Hj in the limit z — > i0 + for i — j and 
fc = I, for any Hamiltonian ([1]). Moreover, for a pure bath 
state, pb — \ip){ip\, the right side of (|3Tj) for i = j and 
k = I can be written as zTi k (z) — Si k z^2 T p i(z) with- 
out any further assumption. The solution of the ensuing 
equation (|33[) is the canonical distribution if the rates Fi k 
satisfy the detailed balance relation. This relation results 
directly from ([3"!?]) as this second-order self-energy leads to 
the Fermi golden rule like expression 

r ik = Y^ e ~ EJT5{E v +ei ~ Ea ~ ek)mu ^ a) \ 2 (38) 

which obviously obeys r k i exp[(efe — £i)/T] — Fi k . The spe- 
cific form ((28)) of the correlations of the bath interaction 
operators U k i is then not decisive to obtain the detailed 
balance relation. 

The crucial result in the derivation of the asymp- 
totic thermal state (l36l) is thus the expression (|30| which 
states that the dynamical correlations (ip\Uij(t)U k i(t')\ip) 
where \ip) is a typical pure bath state of macroscopi- 
cally well-defined energy E, are identical to the respec- 
tive canonical correlation functions with the tempera- 
ture determined by E via (|27p . This result stems from 
(UijUki(T)Ui k (T)Uji)E < exp(iV) for TV » 1 and from 
the equality (|28|) of microcanonical and canonical corre- 
lation functions. These correlation functions are a priori 
related by 

{U 13 U u {t))t = \\ dEn{E)e- E ' T {U l] U k i{T))E (39) 

for any environment. Due to the Boltzmann's relation (|4[), 
the function n(E) exp(—E/T)/Z is essentially a Gaussian 
of variance oc N peaked at the energy E related to T by 
(|27p . For the environment considered in this paper, the 



microcanonical correlations (UijUki(T)) e vary on macro- 
scopic energyscales and hence are equal to canonical av- 
erages as given by (|28)) . 

To better comprehend how the canonical distribution 
for S emerges, we consider the microcanonical form of 
the rate ([55]) . i.e., with Z _1 ^ Q exp(— E a /T) replaced by 
D 1 J2\ a )en E ' ^ e wr ite this expression as 

rik(E) = ± J2 ^ (4°) 

\a)en E 

where the independence of the rate is noted explicitly and 

the definition 7^ = Y,p s i E P - E <* + Lu ik )\{/3\U ik \a)\ 2 is 

used. With these notations and the property llL — U ik , 
we find 

= -p 2^ ilk' = (E) — r ik {E + uj lk ). 

\<x)£H E+u , ih 1 ^ 

(41) 

The last equality follows from D = n(E)8E. It clearly 
shows that the detailed balance relation for the rates Fi k 
results from the Boltzmann's relation (j4|) and from the fact 
that Tifc, which is the Fourier transform of (Ui k U k i(T)) e 
with frequency LOi k , is practically the same for the macro- 
scopic energies E and E + LUi k . 

We finally comment on the weak coupling regime con- 
sidered in the previous section. In this regime, the self- 
energies describing the influence of the bath on S are 
well approximated by their expansions to second order 
in the system-bath interaction Hamiltonian Hi. The de- 
tailed balance relation is satisfied by the rates (|38|) de- 
duced from the second-order self-energies. If we suppose, 
for example, that a part of the Hamiltonian H$ describing 
the intrinsic dynamics of S is of the same order |15j , contri- 
butions stemming from this perturbative self-Hamiltonian 
term must be added to the second-order self-energies. Con- 
sequently, the zeroth-order asymptotic state is no longer 
given by a master equation with rates satisfying the de- 
tailed balance relation. Our results thus indicate, as the 
previous derivations [2, 3 51 4J[5] , that the canonical distribu- 
tion is an attribute of the weak coupling regime. 



6 Conclusion 

In this paper, we have studied the time evolution of a 
system S under the influence of a heat bath of harmonic 
oscillators. The reduced dynamics of S has been derived 
exactly from the Schrodinger equation of the isolated com- 
posite system consisting of S and its bosonic environment. 
Contrary to the usual assumption of a bath initially at 
thermal equilibrium, a pure initial bath state was con- 
sidered. We have shown that, for almost all initial bath 
states of macroscopically well-defined energy, the system 
S relaxes into canonical equilibrium, provided the cou- 
pling to the bath is weak enough. The relaxation of S 
into a canonical mixed state has thus been obtained with- 
out performing any statistical averaging. In other words, 
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no equilibrium state is assumed a priori in our deriva- 
tion. The bath though initially in a pure state induces the 
thermalisation of the system S. The temperature T of the 
asymptotic thermal state of S is determined by the macro- 
scopic energy of the initial bath state. It equals the micro- 
canonical temperature, defined from the bath density of 
states, corresponding to this energy. We have also found 
that, for the overwhelming majority of initial bath states 
of macroscopically well-defined energy, the time evolution 
of S cannot be distinguished from the one obtained assum- 
ing that the boson bath is initially at thermal equilibrium 
at temperature T . 

We emphasize that though these results have been de- 
rived for a boson bath, a substantial part of our approach 
is applicable to other environments. As discussed in the 
previous section, the thermalisation of the system S re- 
sults essentially from two characteristics of the bath: (i) 
the bath density of states verifies the Boltzmann's rela- 
tion (ii) the microcanonical correlation functions of the 
bath operators coupled to S vary on macroscopic ener- 
gyscales. While the first feature is of great generality for 
physical environments, it is not so obvious for the second 
one. It would then be interesting to study the microcanon- 
ical correlations of the coupling operators of other "real- 
istic" environments. Another interesting generalisation of 
the present work would be to add the possibility for the 
system S of interest to exchange also particles with its en- 
vironment and to investigate whether, under these condi- 
tions, a reservoir initially in a "typical" pure state induces 
the relaxation of S into grand-canonical equilibrium. 

We thank R. Chitra for a careful reading of the manuscript. 
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